BID-seq for transcriptome-wide quantitative sequencing of mRNA pseudouridine at base resolution.

Pseudouridine (Ψ) is an abundant RNA modification that is present in and affects the functions of diverse non-coding RNA species, including rRNA, tRNA and small nuclear RNA. Ψ also exists in mammalian mRNA and probably exhibits functional roles; however, functional investigations of mRNA Ψ modifications in mammals have been hampered by the lack of a quantitative method that detects Ψ at base precision. We have recently developed bisulfite-induced deletion sequencing (BID-seq), which provides the community with a quantitative method to map RNA Ψ distribution transcriptome-wide at single-base resolution. Here, we describe an optimized BID-seq protocol for mapping Ψ distribution across cellular mRNAs, which includes fast steps in both library preparation and data analysis. This protocol generates highly reproducible results by inducing high deletion ratios at Ψ modification within diverse sequence contexts, and meanwhile displayed almost zero background deletions at unmodified uridines. When used for transcriptome-wide Ψ profiling in mouse embryonic stem cells, the current protocol uncovered 8,407 Ψ sites from as little as 10 ng of polyA+ RNA input. This optimized BID-seq workflow takes 5 days to complete and includes four main sections: RNA preparation, library construction, next-generation sequencing (NGS) and data analysis. Library construction can be completed by researchers who have basic knowledge and skills in molecular biology and genetics. In addition to the experimental protocol, we provide BID-pipe ( https://github.com/y9c/pseudoU-BIDseq ), a user-friendly data analysis pipeline for Ψ site detection and modification stoichiometry quantification, requiring only basic bioinformatic and computational skills to uncover Ψ signatures from BID-seq data.

Transcriptome-wide analysis of pseudouridylation in Drosophila melanogaster.

Pseudouridine (Psi) is one of the most frequent post-transcriptional modification of RNA. Enzymatic Psi modification occurs on rRNA, snRNA, snoRNA, tRNA, and non-coding RNA and has recently been discovered on mRNA. Transcriptome-wide detection of Psi (Psi-seq) has yet to be performed for the widely studied model organism Drosophila melanogaster. Here, we optimized Psi-seq analysis for this species and have identified thousands of Psi modifications throughout the female fly head transcriptome. We find that Psi is widespread on both cellular and mitochondrial rRNAs. In addition, more than a thousand Psi sites were found on mRNAs. When pseudouridylated, mRNAs frequently had many Psi sites. Many mRNA Psi sites are present in genes encoding for ribosomal proteins, and many are found in mitochondrial encoded RNAs, further implicating the importance of pseudouridylation for ribosome and mitochondrial function. The 7SLRNA of the signal recognition particle is the non-coding RNA most enriched for Psi. The 3 mRNAs most enriched for Psi encode highly expressed yolk proteins (Yp1, Yp2, and Yp3). By comparing the pseudouridine profiles in the RluA-2 mutant and the w1118 control genotype, we identified Psi sites that were missing in the mutant RNA as potential RluA-2 targets. Finally, differential gene expression analysis of the mutant transcriptome indicates a major impact of loss of RluA-2 on the ribosome and translational machinery.

Pseudouridines in RNAs: switching atoms means shifting paradigms.

The structure, stability, and function of various coding and noncoding RNAs are influenced by chemical modifications. Pseudouridine (Ψ) is one of the most abundant post-transcriptional RNA base modifications and has been detected at individual positions in tRNAs, rRNAs, mRNAs, and snRNAs, which are referred to as Ψ-sites. By allowing formation of additional bonds with neighboring atoms, Ψ strengthens RNA-RNA and RNA-protein interactions. Although many aspects of the underlying modification reactions remain unclear, the advent of new transcriptome-wide methods to quantitatively detect Ψ-sites has recently changed our perception of the functional roles and importance of Ψ. For instance, it is now clear that the occurrence of Ψs appears to be directly linked to the lifetime and the translation efficiency of a given mRNA molecule. Furthermore, the administration of Ψ-containing RNAs reduces innate immune responses, which appears strikingly advantageous for the development of generations of mRNA-based vaccines. In this review, we aim to comprehensively summarize recent discoveries that highlight the impact of Ψ on various types of RNAs and outline possible novel biomedical applications of Ψ.

Pseudouridine as a novel biomarker in prostate cancer.

Epitranscriptomic analysis has recently led to the profiling of modified nucleosides in cancer cell biological matrices, helping to elucidate their functional roles in cancer and reigniting interest in exploring their use as potential markers of cancer development and progression. Pseudouridine, one of the most well-known and the most abundant of the RNA nucleotide modifications, is the C5-glycoside isomer of uridine and its distinctive physiochemical properties allows it to perform many essential functions. Pseudouridine functionally (a) confers rigidity to local RNA structure by enhancing RNA stacking, engaging in a cooperative effect on neighboring nucleosides that overall contributes to RNA stabilization (b) refines the structure of tRNAs, which influences their decoding activity (c) facilitates the accuracy of decoding and proofreading during translation and efficiency of peptide bond formation, thus collectively improving the fidelity of protein biosynthesis and (e) dynamically regulates mRNA coding and translation. Biochemical synthesis of pseudouridine is carried out by pseudouridine synthases. In this review we discuss the evidence supporting an association between elevated pseudouridine levels with the incidence and progression of human prostate cancer and the translational significance of the value of this modified nucleotide as a novel biomarker in prostate cancer progression to advanced disease.

Direct Determination of Pseudouridine in RNA by Mass Spectrometry Coupled with Stable Isotope Labeling.

Pseudouridine (Ψ) is the only "mass-silent" nucleoside produced by post-transcriptional RNA modification. We developed a mass spectrometry (MS)-based technique coupled with in vivo deuterium (D) labeling of uridines for direct determination of Ψs in cellular RNA and applied it to the comprehensive analysis of post-transcriptional modifications in human ribosomal RNAs. The method utilizes human TK6/mouse FM3A cells deficient in uridine monophosphate synthase using a CRISPR-Cas9 technique to turn off de novo uridine synthesis and fully labels uridines with D at uracil positions 5 and 6 by cultivating the cells in a medium containing uridine-5,6-D2. The pseudouridylation reaction in those cells results in the exchange of the D at the C5 of uracil with hydrogen from solvent, which produces a -1 Da mass shift, thus allowing MS-based determination of RNA Ψs. We present here the experimental details of this method and show that it allows the identification of all Ψs in human major nuclear and nucleolar RNAs, including several previously unknown Ψs. Because the method allows direct determination of Ψs at the femtomole level of RNA, it will serve as a useful tool for structure/function studies of a wide variety of noncoding RNAs.

Potential G-Quadruplex Forming Sequences and N6-Methyladenosine Colocalize at Human Pre-mRNA Intron Splice Sites.

Maturation of mRNA in humans involves modifying the 5' and 3' ends, splicing introns, and installing epitranscriptomic modifications that are essential for mRNA biogenesis. With respect to epitranscriptomic modifications, they are usually installed in specific consensus motifs, although not all sequences are modified suggesting a secondary structural component to site selection. Using bioinformatic analysis of published data, we identify in human mature-mRNA that potential RNA G-quadruplex (rG4) sequences colocalize with the epitranscriptomic modifications N6-methyladenosine (m6A), pseudouridine (Ψ), and inosine (I). Using the only available pre-mRNA data sets from the literature, we demonstrate colocalization of potential rG4s and m6A was greatest overall and occurred in introns near 5' and 3' splice sites. The loop lengths and sequence context of the m6A-bearing potential rG4s exhibited short loops most commonly comprised of single A nucleotides. This observation is consistent with a literature report of intronic m6A found in SAG (S = C or G) consensus motifs that are also recognized by splicing factors. The localization of m6A and potential rG4s in pre-mRNA at intron splice junctions suggests that these features could function together in alternative splicing. A similar analysis for potential rG4s around sites of Ψ installation or A-to-I editing in mRNA also found a colocalization; however, the frequency was less than that observed with m6A. These bioinformatic analyses guide a discussion of future experiments to understand how noncanonical rG4 structures may collaborate with epitranscriptomic modifications in the human cellular context to impact cellular phenotype.

RNA modifications and the link to human disease.

Ribonucleic acid (RNA) is involved in translation and transcription, which are the mechanisms in which cells express genes (Alberts et al., 2002). The three classes of RNA discussed are transfer RNA (tRNA), messenger RNA (mRNA), and ribosomal RNA (rRNA). mRNA is the transcript encoded from DNA, rRNA is associated with ribosomes, and tRNA is associated with amino acids and is used to read mRNA transcripts to make proteins (Lodish, Berk, Zipursky, et al., 2000). Interestingly, the function of tRNA, rRNA, and mRNA can be significantly altered by chemical modifications at the co-transcriptional and post-transcriptional levels, and there are over 171 of these modifications identified thus far (Boccaletto et al., 2018; Modomics-Modified bases, 2017). Several of these modifications are linked to diseases such as cancer, diabetes, and neurological disorders. In this review, we will introduce a few RNA modifications with biological functions and how dysregulation of these RNA modifications is linked to human disease.

A sensitive and radiolabeling-free method for pseudouridine detection.

Existing methodologies for detecting Pseudouridine (Ψ) mostly use CMCT labeling or radiolabeling. Described herein is a sensitive and quantitative method for Ψ detection that does not need this labelling. This approach combines the selectivity of a 10-23 DNAzyme, which can distinguish Ψ from uridine (U), with rolling circle amplification (RCA) to increase the sensitivity of the assay.

PseUI: Pseudouridine sites identification based on RNA sequence information.

BACKGROUND: Pseudouridylation is the most prevalent type of posttranscriptional modification in various stable RNAs of all organisms, which significantly affects many cellular processes that are regulated by RNA. Thus, accurate identification of pseudouridine (Ψ) sites in RNA will be of great benefit for understanding these cellular processes. Due to the low efficiency and high cost of current available experimental methods, it is highly desirable to develop computational methods for accurately and efficiently detecting Ψ sites in RNA sequences. However, the predictive accuracy of existing computational methods is not satisfactory and still needs improvement. RESULTS: In this study, we developed a new model, PseUI, for Ψ sites identification in three species, which are H. sapiens, S. cerevisiae, and M. musculus. Firstly, five different kinds of features including nucleotide composition (NC), dinucleotide composition (DC), pseudo dinucleotide composition (pseDNC), position-specific nucleotide propensity (PSNP), and position-specific dinucleotide propensity (PSDP) were generated based on RNA segments. Then, a sequential forward feature selection strategy was used to gain an effective feature subset with a compact representation but discriminative prediction power. Based on the selected feature subsets, we built our model by using a support vector machine (SVM). Finally, the generalization of our model was validated by both the jackknife test and independent validation tests on the benchmark datasets. The experimental results showed that our model is more accurate and stable than the previously published models. We have also provided a user-friendly web server for our model at http://zhulab.ahu.edu.cn/PseUI , and a brief instruction for the web server is provided in this paper. By using this instruction, the academic users can conveniently get their desired results without complicated calculations. CONCLUSION: In this study, we proposed a new predictor, PseUI, to detect Ψ sites in RNA sequences. It is shown that our model outperformed the existing state-of-art models. It is expected that our model, PseUI, will become a useful tool for accurate identification of RNA Ψ sites.

Transcription of Leishmania major U2 small nuclear RNA gene is directed by extragenic sequences located within a tRNA-like and a tRNA-Ala gene.

BACKGROUND: Leishmania and other trypanosomatid parasites possess atypical mechanisms of gene expression, including the maturation of mRNAs by trans-splicing and the involvement of RNA Polymerase III in transcription of all snRNA molecules. Since snRNAs are essential for trans-splicing, we are interested in the study of the sequences that direct their expression. Here we report the characterization of L. major U2 snRNA promoter region. RESULTS: All species of Leishmania possess a single U2 snRNA gene that contains a divergently-oriented tRNA-Ala gene in the upstream region. Between these two genes we found a tRNA-like sequence that possesses conserved boxes A and B. Primer extension and RT-qPCR analyses with RNA from transiently-transfected cells showed that transcription of L. major U2 snRNA is almost abolished when boxes A and B from the tRNA-like are deleted or mutated. The levels of the U2 snRNA were also highly affected when base substitutions were introduced into box B from the tRNA-Ala gene and the first nucleotides of the U2 snRNA gene itself. We also demonstrate that the tRNA-like is transcribed, generating a main transcript of around 109 bases. As pseudouridines in snRNAs are required for splicing in other organisms, we searched for this modified nucleotide in the L. major U2 snRNA. Our results show the presence of six pseudouridines in the U2 snRNA, including one in the Sm site that has not been reported in other organisms. CONCLUSIONS: Four different regions control the transcription of the U2 snRNA gene in L. major: boxes A and B from the neighbor tRNA-like, box B from the upstream tRNA-Ala gene and the first nucleotides of the U2 snRNA. Thus, the promoter region of L. major U2 snRNA is different from any other promoter reported for snRNAs. Pseudouridines could play important roles in L. major U2 snRNA, since they were found in functionally important regions, including the branch point recognition region and the Sm binding site.

Transcriptome-Wide Identification of Pseudouridine Modifications Using Pseudo-seq.

A diverse array of post-transcriptional modifications is found in RNA molecules from all domains of life. While the locations of RNA modifications are well characterized in abundant noncoding RNAs, modified sites in less abundant mRNAs are just beginning to be discovered. Recent work has revealed hundreds of previously unknown and dynamically regulated pseudouridines (Ψ) in mRNAs from diverse organisms. This unit describes Pseudo-seq, an efficient, high-resolution method for identification of Ψs genome-wide. This unit includes methods for isolation of RNA from S. cerevisiae, preparation of Pseudo-seq libraries from RNA samples, and identification of sites of pseudouridylation from the sequencing data. Pseudo-seq is applicable to any organism or cell type, facilitating rapid identification of novel pseudouridylation events.

PPUS: a web server to predict PUS-specific pseudouridine sites.

MOTIVATION: Pseudouridine (Ψ), catalyzed by pseudouridine synthase (PUS), is the most abundant RNA modification and has important cellular functions. Developing an algorithm to identify Ψ sites is an important work. And it is better if the algorithm could assign which PUS modifies the Ψ sites. Here, we developed PPUS (http://lyh.pkmu.cn/ppus/), the first web server to predict PUS-specific Ψ sites. PPUS: employed support vector machine as the classifier and used nucleotides around Ψ sites as the features. Currently, PPUS: could accurately predict new Ψ sites for PUS1, PUS4 and PUS7 in yeast and PUS4 in human. PPUS: is well designed and friendly to user. AVAILABILITY AND IMPLEMENTATION: Our web server is available freely for non-commercial purposes at: http://lyh.pkmu.cn/ppus/ CONTACT: liyanhui@bjmu.edu.cn or cuiqinghua@hsc.pku.edu.cn.

Pseudouridine in a new era of RNA modifications.

Two articles recently published in Nature and Cell report the first transcriptome-wide maps of pseudouridine (Ψ) at single-base resolution through selective chemical labeling, suggesting new mechanisms and functions of Ψ in mRNA and non-coding RNA molecules.

Transcriptome-wide mapping of pseudouridines: pseudouridine synthases modify specific mRNAs in S. cerevisiae.

We developed a novel technique, called pseudouridine site identification sequencing (PSI-seq), for the transcriptome-wide mapping of pseudouridylation sites with single-base resolution from cellular RNAs based on the induced termination of reverse transcription specifically at pseudouridines following CMCT treatment. PSI-seq analysis of RNA samples from S. cerevisiae correctly detected all of the 43 known pseudouridines in yeast 18S and 25S ribosomal RNA with high specificity. Moreover, application of PSI-seq to the yeast transcriptome revealed the presence of site-specific pseudouridylation within dozens of mRNAs, including RPL11a, TEF1, and other genes implicated in translation. To identify the mechanisms responsible for mRNA pseudouridylation, we genetically deleted candidate pseudouridine synthase (Pus) enzymes and reconstituted their activities in vitro. These experiments demonstrated that the Pus1 enzyme was necessary and sufficient for pseudouridylation of RPL11a mRNA, whereas Pus4 modified TEF1 mRNA, and Pus6 pseudouridylated KAR2 mRNA. Finally, we determined that modification of RPL11a at Ψ -68 was observed in RNA from the related yeast S. mikitae, and Ψ -239 in TEF1 mRNA was maintained in S. mikitae as well as S. pombe, indicating that these pseudouridylations are ancient, evolutionarily conserved RNA modifications. This work establishes that site-specific pseudouridylation of eukaryotic mRNAs is a genetically programmed RNA modification that naturally occurs in multiple yeast transcripts via distinct mechanisms, suggesting that mRNA pseudouridylation may provide an important novel regulatory function. The approach and strategies that we report here should be generally applicable to the discovery of pseudouridylation, or other RNA modifications, in diverse biological contexts.

Pseudouridine profiling reveals regulated mRNA pseudouridylation in yeast and human cells.

Post-transcriptional modification of RNA nucleosides occurs in all living organisms. Pseudouridine, the most abundant modified nucleoside in non-coding RNAs, enhances the function of transfer RNA and ribosomal RNA by stabilizing the RNA structure. Messenger RNAs were not known to contain pseudouridine, but artificial pseudouridylation dramatically affects mRNA function--it changes the genetic code by facilitating non-canonical base pairing in the ribosome decoding centre. However, without evidence of naturally occurring mRNA pseudouridylation, its physiological relevance was unclear. Here we present a comprehensive analysis of pseudouridylation in Saccharomyces cerevisiae and human RNAs using Pseudo-seq, a genome-wide, single-nucleotide-resolution method for pseudouridine identification. Pseudo-seq accurately identifies known modification sites as well as many novel sites in non-coding RNAs, and reveals hundreds of pseudouridylated sites in mRNAs. Genetic analysis allowed us to assign most of the new modification sites to one of seven conserved pseudouridine synthases, Pus1-4, 6, 7 and 9. Notably, the majority of pseudouridines in mRNA are regulated in response to environmental signals, such as nutrient deprivation in yeast and serum starvation in human cells. These results suggest a mechanism for the rapid and regulated rewiring of the genetic code through inducible mRNA modifications. Our findings reveal unanticipated roles for pseudouridylation and provide a resource for identifying the targets of pseudouridine synthases implicated in human disease.

Transcriptome-wide mapping reveals widespread dynamic-regulated pseudouridylation of ncRNA and mRNA.

Pseudouridine is the most abundant RNA modification, yet except for a few well-studied cases, little is known about the modified positions and their function(s). Here, we develop Ψ-seq for transcriptome-wide quantitative mapping of pseudouridine. We validate Ψ-seq with spike-ins and de novo identification of previously reported positions and discover hundreds of unique sites in human and yeast mRNAs and snoRNAs. Perturbing pseudouridine synthases (PUS) uncovers which pseudouridine synthase modifies each site and their target sequence features. mRNA pseudouridinylation depends on both site-specific and snoRNA-guided pseudouridine synthases. Upon heat shock in yeast, Pus7p-mediated pseudouridylation is induced at >200 sites, and PUS7 deletion decreases the levels of otherwise pseudouridylated mRNA, suggesting a role in enhancing transcript stability. rRNA pseudouridine stoichiometries are conserved but reduced in cells from dyskeratosis congenita patients, where the PUS DKC1 is mutated. Our work identifies an enhanced, transcriptome-wide scope for pseudouridine and methods to dissect its underlying mechanisms and function.

Absolute and relative quantification of RNA modifications via biosynthetic isotopomers.

In the resurging field of RNA modifications, quantification is a bottleneck blocking many exciting avenues. With currently over 150 known nucleoside alterations, detection and quantification methods must encompass multiple modifications for a comprehensive profile. LC-MS/MS approaches offer a perspective for comprehensive parallel quantification of all the various modifications found in total RNA of a given organism. By feeding (13)C-glucose as sole carbon source, we have generated a stable isotope-labeled internal standard (SIL-IS) for bacterial RNA, which facilitates relative comparison of all modifications. While conventional SIL-IS approaches require the chemical synthesis of single modifications in weighable quantities, this SIL-IS consists of a nucleoside mixture covering all detectable RNA modifications of Escherichia coli, yet in small and initially unknown quantities. For absolute in addition to relative quantification, those quantities were determined by a combination of external calibration and sample spiking of the biosynthetic SIL-IS. For each nucleoside, we thus obtained a very robust relative response factor, which permits direct conversion of the MS signal to absolute amounts of substance. The application of the validated SIL-IS allowed highly precise quantification with standard deviations<2% during a 12-week period, and a linear dynamic range that was extended by two orders of magnitude.

Mass spectrometry-based quantification of pseudouridine in RNA.

Direct detection of pseudouridine (ψ), an isomer of uridine, in RNA is challenging. The most popular method requires chemical derivatization using N-cyclohexyl-N'-ß-(4-methylmorpholinum ethyl) carbodiimide p-tosylate (CMCT) followed by radiolabeled primer extension mediated by reverse transcriptase. More recently, mass spectrometry (MS)-based approaches for sequence placement of pseudouridine in RNA have been developed. Nearly all of these approaches, however, only yield qualitative information regarding the presence or absence of pseudouridine in a given RNA population. Here, we have extended a previously developed liquid chromatography tandem mass spectrometry (LC-MS/MS) method to enable both the qualitative and quantitative analysis of pseudouridine. Quantitative selected reaction monitoring (SRM) assays were developed using synthetic oligonucleotides, with or without pseudouridine, and the results yielded a linear relationship between the ion abundance of the pseudouridine-specific fragment ion and the amount of pseudouridine-containing oligonucleotide present in the original sample. Using this quantitative SRM assay, the extent of pseudouridine hypomodification in the conserved T-loop of tRNA isolated from two different Escherichia coli strains was established.

The ribosome assembly factor Nep1 responsible for Bowen-Conradi syndrome is a pseudouridine-N1-specific methyltransferase.

Nep1 (Emg1) is a highly conserved nucleolar protein with an essential function in ribosome biogenesis. A mutation in the human Nep1 homolog causes Bowen-Conradi syndrome-a severe developmental disorder. Structures of Nep1 revealed a dimer with a fold similar to the SPOUT-class of RNA-methyltransferases suggesting that Nep1 acts as a methyltransferase in ribosome biogenesis. The target for this putative methyltransferase activity has not been identified yet. We characterized the RNA-binding specificity of Methanocaldococcus jannaschii Nep1 by fluorescence- and NMR-spectroscopy as well as by yeast three-hybrid screening. Nep1 binds with high affinity to short RNA oligonucleotides corresponding to nt 910-921 of M. jannaschii 16S rRNA through a highly conserved basic surface cleft along the dimer interface. Nep1 only methylates RNAs containing a pseudouridine at a position corresponding to a previously identified hypermodified N1-methyl-N3-(3-amino-3-carboxypropyl) pseudouridine (m1acp3-Psi) in eukaryotic 18S rRNAs. Analysis of the methylated nucleoside by MALDI-mass spectrometry, HPLC and NMR shows that the methyl group is transferred to the N1 of the pseudouridine. Thus, Nep1 is the first identified example of an N1-specific pseudouridine methyltransferase. This enzymatic activity is also conserved in human Nep1 suggesting that Nep1 is the methyltransferase in the biosynthesis of m1acp3-Psi in eukaryotic 18S rRNAs.